The structure of intercalated water in superconducting Nao.35Co02 l.37D20: 
Implications for the superconducting phase diagram 
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We have used electron and neutron powder diffraction to elucidate the structural properties of 
superconducting Naa:Co02 -2/020. Our measurements show that our superconducting sample exhbits 
a number of supercells ranging from ia* to j^a*, but the most predominant one, observed also in 
the neutron data, is a double hexagonal cell with dimensions 2a x 2a x c. Rietveld analysis reveals 
that D2O is inserted between C0O2 sheets as to form a layered network of NaOe triangular prisms. 
Our model removes the need to invoke a 5K superconducting point compound and suggests that a 
solid solution of Na is possible within a constant amount of water y. 



The tuning of Na content in the alkali layered cobal- 
tate Na2;Co02 results in remarkable changes in its physi- 
cal behavior, ranging from magneto-thermoelectricity at 
a;=O.750| to charge ordering at a;=O.5[3,0 and 5K super- 
conductivity at x=0.3 after intercalation of H2O. 4J This 
interesting combination of physical phenomena occurs 
within a structural motif reminiscent of geometrically 
frustrated systems. Here, C0O2 sheets are constructed 
by edge sharing CoOe octahedra, forming a quasi-2D tri- 
angular net. The exact mechanism of superconductiv- 
ity within this triangular motif is currently the topic of 
extensive scientific inquiry, as it holds the possibility of 
realizing Anderson's Resonating Valence Bond model. 

Water insertion increases the separation between C0O2 
sheets and band structure calculations suggest that 
this effectively minimizes the electronic inter-planar 
coupling. H 0This dimensional cross-over effect is cor- 
roborated by measurements showing that Tn increases 
with the separation between C0O2 sheets. However, 
the same measurements also show that Tc does not vary 
significantly with electronic doping x. This contradicts 
a proposed analogy between this alkali hydrate and the 
high-Tc cuprates [lOj , whereby optimal Tc would occur 
within a narrow range of electronic doping. 

The resolution of this controversy may lie in the details 
of the crystal structure of this alkali hydrate. Two mutu- 
ally incompatible structural models of water-intercalated 
of Naa;Co02 • yH20 , have been proposed to date.flllll^ 
The model proposed by Lynn et al. asserts that the inter- 
calated water assumes a structure similar to zce^3| in be- 
tween the C0O2 sheets, whereas the one by Jorgensen et 
al. suggests a coordination and crystal chemical linkage 
between Na and H2O similar that found in many alkali 
hydrate systems. 12] Both structural models imply that 
Na concentration x and the amount of intercalated wa- 
ter y are intimately related. But Jorgensen's model also 
implies that optimal Tc corresponds to a line compound 



exhibiting full Na and D2O ordering with a Na/H2 ratio 
~l/4. J^. Both models are based on the parent structure 
and unit cell, and rely heavily on chemical constraints to 
interpret the multiple symmetry-equivalent water sites 
in that cell. In spite of their complexity, these models 
do not reproduce some important aspects of the data. 
In particular, an intense Bragg peak at (2.8A), and a 
broad distribution of intensity around (^2.6A) are not 
accounted for, although Jorgensen et al. speculate that 
both features may result from a larger periodicity. 

To unravel the complexity of its crystal structure and 
gain insight into parameters that are critical in estab- 
lishing a superconducting phase diagram, we clearly 
need to establish the true translational symmetry of 
Naa;Co02 • J/H2O . To this end, we have performed elec- 
tron diffraction (ED) and neutron powder diffraction 
(NPD) measurements, the latter employing isotopic H- 
D substitutions. Both sets of data clearly establish that 
the dominant ordering mode is a doubling of the hexag- 
onal unit cell along both a and b. Our analysis shows 
that Na-atoms coordinate with D2O to form an ordered 
network of NaOg triangular prisms. This model can ac- 
count well for the intense Bragg reflection in the neutron 
data at ~2.8A. ED also shows a number of different su- 
perstructures, suggesting significant Na-inhomogeneities. 
These additional ordering modes are most likely the ori- 
gin the of a broad distribution of scattered intensity at 
^2.6 A, a common feature of all the neutron data pub- 
lished thus far. Our double-cell model is able to accom- 
modate a variable amount of Na within the same water 
framework, strongly suggesting that Naa;Co02 • yH20 is 
not a line compound. Our evidence of more than one 
Na/water ordering mode within a presumably coherent 
C0O2 network should open a serious debate about the 
true nature of the superconducting phase. 

Description of the synthesis and characterization of 
these samples is found elsewhere. |1| In this work two sam- 
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pies were investigated with a;=0.35, which exhibited a su- 
perconducting transitions at 4.5K and 2.8K.0| ED pat- 
terns of these samples were taken using a Philips CM30 
transmission electron microscope operated at 300kV with 
the sample cooled to approximately 90K using a Gatan 
double tilt liquid nitrogen stage. NPD data were mea- 
sured from these polycrystalline samples at 2K using the 
high resolution powder diffractometer E9 (A=1.7973A), 
located at the Berlin Neutron Scattering Center, at the 
Hahn-Meitner-Institut. Additional data were collected 
on samples of varying x on the GEM diffractometer at 
ISIS which were critical in developing the present model, 
and will be presented elsewhere. To located the 0-atoms 
of the water molecules, the sample with Tc=2.8K that 
was originally hydrated with D2O, was exposed to a 
(H20)o.64i (D20)o.359 moisturc mixture at a tempera- 
ture of 26°C for 3 days. This isotopic H/D ratio gives a 
zero mean coherent neutron scattering length for protons 
(&//=-0.374,6d=0.667), so in principle for this sample the 
contribution of protons to the entire diffraction pattern 
can be ignored. 

Our ED measurements show that the average symme- 
try of the superconducting compound remains hexagonal 
everywhere. In the majority of grains, super-lattice re- 
flections indicate a cell doubling giving an average unit 
cell of 2a X 2a X c (see fig. ^) where a and c are the 
lattice constants of the parent VQ^/mmc hexagonal cell. 
Examination of the NPD data shows the strong reflection 
observed at 37°(2.8A) can be assigned to the (|,^,0), 
(i, i, 1) Bragg peaks, as shown in fig. We can con- 
struct a chemically sensible double-cell structural model 
that fits the 2.8A peak, providing the strongest indica- 
tion that this and the ED superlattice reflections have 
the same origin. Further investigation of the ED data 
reveals an even more complex superstructure in some ar- 
eas of the sample as seen in fig. Satellite reflections, 
decorating the parent hexagonal lattice reflections, are 
arranged in hexagonal nets with q- vectors whose magni- 
tude has been seen to vary from l/15a* to l/3a*, but 
whose direction is always parallel or nearly parallel to 
< 110 >*. The example shown in fig corresponds 
to a supercell translation vector of 3a + 3|6, requiring 
a discommensuration to allow registry between the su- 
percells, equivalent to the Nal-Na2 interatomic vector 
(see below). In all such patterns mirror symmetry per- 
pendicular to the basal planes is broken as seen by the 
gross asymmetry in the intensities of satellite reflections 
of type g ± q. However, the 6-fold symmetry appears to 
be conserved throughout. The detailed superstructure 
of these modulated phases will be discussed further in a 
forthcoming publication. 

These more complex super-structures seen in the ED 
measurements are not clearly observed in our NPD data 
Indeed, an examination of the NPD data does not show 
the presence of any additional super-lattice reflections 
above 37° suggesting that disorder (which we model as a 
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FIG. 1: (a-b) ED patterns from Nao.35CoO2.yD2O. The 
arrows in (a) indicate the cell doubling reflections. (c) 
Neutron powder diffraction data measured at 2K from 
Nao.35Co02-yD20. The difference between the observations 
and the model are plotted under the diffraction pattern, 
while the expected location of reflection are shown as ver- 
tical lines in between. To account for the broad scatter- 
ing at centered at 40° (indicated by arrow in inset) the 
background was interpolated manually for the entire pat- 
tern. The inset in (c) depicts a portion of the neutron 
diffractogram from the Nao.35Co02 • J/D2O (heavy line) and 
Nao.35Co02 • y(Ho.64iDo. 359)20 (light line). The background 
arising from the incoherent neutron scattering from H-atoms 
has been subtracted. 



large effective Debye- Waller factor for water) is substan- 
tially attenuating their intensity (superstructure reflec- 
tions in ED are enhanced due to dynamical diffraction). 
Comparison of the diffractograms from our deuterated 
and isotopically mixed HD samples show appreciable dif- 
ferences as illustrated in the inset of fig. While some 
reflections show no significant change in their intensity 
with the isotopic substitution (for example 29 ~31 and 
37°), others do. That the intensity of the super-lattice 
reflections at 37° as well as the broad feature at 40° is 
changed, indicates that their structure factor has an ap- 
preciable contribution from protons. The origin of the 
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TABLE I: Structural parameters determined from Rietveld anal- 
ysis of our NPD data measured from our deuterated x = 0.35 
sample at 2K. The refined lattice constants were a =5. 63477(21)A 
and c =19.5456(13)A (space group PGs/m). The weighted profile 
R-factor for this refinement is •wRp=fi.22%. The occupancies of the 
Na sites were constrained so that f{Nal) + f(Na2) = 0.35. The 
deuterons labeled D31 and D32 are coordinated to 03, and simi- 
larly for 04,D41,D42. The refined total amount of D2O was com- 
puted to be j;=1.37(l). In the Rietveld analysis D2O was inserted 
as a rigid body, with a D-0 bond-length constrained to 0.99 A and 
the D-O-D torsion angle set to 108°. C/iso(Co,O,Na)=0.0005(2)A2 
and t/iso(D2O)=0.035{4)A2. Selected bondlengths (in A) are given 
at the bottom of the table. 
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broad diffuse feature at 40° (2.6A) is controversial as it 
has been described eitlier as an impurity phase [Tlj or a 
broad Bragg reflectionfl^. We argue below that it may 
arise from a superposition of Bragg peaks from the addi- 
tional ordering modes and from short-range order. 

To construct a structural model of the Naj;Co02-yD2 
superconductor we analyzed the NPD data using a dou- 
ble hexagonal unit cell (2a x 2a x c). The space group 
was identified as PBs/m using symmetry analysis of the 
parent PBs/mmc space group and taking into consider- 
ation the ED observations. The solution of the crystal 
structure was achieved by first analyzing the isotopically 
exchanged H-D sample to located the positions of the O- 
atoms of the water molecules, and then the deuterated 
sample to located the positions of the protons. This ap- 
proach produced a good agreement between our model 
and the NPD data, including the intensity of the reflec- 
tions at 37°as shown in fig. Structural parameters 
and occupancies determined from the Rietveld analysis 
of the deuterated sample are given, in table 

The main structural features of our model of supercon- 
ducting Naj;Co02 • yD20 are shown in figures |21 and 13 
The structure of C0O2 sheets does not differ substantially 




FIG. 2: (Color) A three dimensional representation of the struc- 
ture of Na2,Co02 ■ iyD20 obtained from our Rietveld analysis. For 
simplicity we show only half the unit cell along the c-axis and only 
the Nal site with the D2O molecule centered at 03. The Na-03 
and D-O(-Co) bond-lengths are shown. 

from the one proposed by Takada. Q . To simplify the dis- 
cussion of the Na and D2O environment, we first consider 
an ideal structure of composition x=^ and y=1.5 con- 
sisting of a fully occupied Nal site surrounded by 6D2O 
molecules in a triangular prism geometry (see fig. |3^). 
Each NaOg prism consists of a Na-atom in the center 
of the prism and 0-atoms at each corner above and be- 
low the Na with Na-0 bond-lengths of 2.4A (fig. El- For 
this composition, NaOg prisms would be uncoupled from 
each other. This is a very typical coordination for Na 
in Na- hydrate systems. [13 In such materials Na atoms 
are surrounded by H2O forming hydration shells, with 
the 0-atom bonded to the alkali metal with a typical 
bond length of 2.4A, while the protons are pointing away 
from it.[13| The most common Na-H20 coordination is 
6, close to the value we find here. The doubling of a and 
b originates from this ordered arrangement of Na and its 
hydration shell (see fig. OJ- In fact, our NPD analysis 
has identified two Na sites centered at Nal (0,0,^) and 
Na2 (i,|,i). It is easy to see that additional Na (up to 
x—^) with the same coordination number can be accom- 
modated into the ideal structure by the progressive filling 
of the Na2 sites, without a substantial rearrangement of 
the water molecules. We find that the Nal site is almost 
fully occupied with fractional occupancy / ~ 0.88(4), 
while the second Na2 site is significantly less populated 
with / ^0.56(4). The coordination of 0-atoms around 
the Na2 site is naturally the same as in the Nal, such 
that NaOg prisms can be centered at either Na site. The 
total amount of D2O determined from our Rietveld re- 
finements is y =1.37(1), in good agreement with other 
reports. |j, |lj |l^ |lj| Finally our analysis shows that at 
least one of the protons of the D2O molecule tend to 
point toward the C0O2 sheets forming a hydrogen bond 
(D-0=1.7 A) with an 0-atom in that layer (fig.Ell. The 
structural model we employed in our best fits is slightly 
more complex than the one we just described, in that 



it includes a small amount of additional water centered 
around the 04 position (see table). This position, equiva- 
lent to the 03 position in the parent unit cell, may repre- 
sent contributions of the four possible coherent domains 
of the double hexagonal cell. More likely, however, this 
site helps modeling the large amount of disorder in the 
structure (see below). 

The main implication of our model, and the most rele- 
vant for the physics of this compound, is that Naa;Co02 • 
2/H2O has the potential of being a homogeneous solid so- 
lution in a wide range of Na concentrations, rather than a 
point compound as suggested by Jorgensen et al. 12]. A 
second important implication is that changes of Na con- 
tent can, within limits, be accommodated at constant 
water content, as each 0-atom of the D2O molecule lo- 
cally will be coordinated two Na-atoms. In reality, a 
series of more or less stable ordered arrangements of Na2 
occupied and vacant sites will occur for specific Na con- 
centrations, although not all of them will be kinetically 
accessible. The complex synthesis route of this mate- 
rial (de-intercalation of Na and intercalation of water) 
helps to explain why more than one superstructure type 
is observed in the same sample. The inhomogeneous dis- 
tribution of Na that is expected from de-intercalation 
produces a disordered background on which intercalated 
H2O must arrange. The intercalation of water would 
yield some re-organization of the Na-atoms at the local 
scale, but energetically a global re-arrangement is too 
costly (unless an anneal is possible) . As water is inserted 
into the lattice, local ordering of Na atoms occurs that 
resembles as much as possible the structure we describe. 
Fluctuations in the Na concentration are accommodated 
locally into the Na2-site, occasionally reaching an ap- 
propriate value for the formation of a long-period super- 
structure. This model would account for the predomi- 
nant doubling mode and the many other diverse super- 
structures observed by ED. We note that our model does 
fold onto the structural models proposed thus far. flllll^ 
For example the Nal site of the model by Lynn et al. 
maps onto the Nal site of our model and similarly for 
the model of Jorgensen et al. and the Na2 site. 

The presence of these inhomogeneities should lead one 
to reconsider carefully the true nature of the supercon- 
ducting phase. For example the predominant 2a x 2a x c 
ordering may represent the 5K superconducting point 
phase while other supercells may also be superconduct- 
ing with a lower Tc. This scenario is consistent with the 
almost constant variation of Tc that has been observed 
as function x-J^ NPD data over a range of x that we have 
measured (0.28< x <0.37) display the (i, 5,0)-(i, i, 1) 
super-lattice reflection pair, indicating that the predom- 
inate of this ordering is quite stable over this x-range. 
Other scenarios are possible, such as those in which a 
minority phase would be superconducting, and would 
partially shield the majority non-superconducting phase, 
but this is unlikely as superconducting fractions are high. 
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FIG. 3: (Color) Projections of the Na-D20 layer down the c— axis. 
The NaOe triangular prisms centered on Na are highlighted. The 
2a X 2a X c unit cell is shown in solid lines and the a X a X c cell 
as dashed lines. The doubling of the unit cell is evident from the 
position of the Na sites. D-atoms are omitted for clarity, (a) The 
idealized x = j structure, (b) A model for the a; = i structure 
with the Na2 site filled. The coordination to the 04 site is shown. 



indicating that a majority phase is a superconductor. 

In summary, our neutron and electron diffraction data 
have unraveled a new aspect of the structure of Naa;Co02- 
J/D2O , shedding light into its material and physical prop- 
erties. The formation of a NaOg triangular prism net- 
work can sustain a large range of x, without the neces- 
sity to vary the water content, or the separation between 
C0O2 sheets. We argue on the basis of the structure that 
the amount of intercalated water essentially defines the 
separation between C0O2 sheets and thus Tc. A variety 
of sodium dopings, some of which leads to additional or- 
dering can be accommodated in the structure and should 
have little influence on separation of C0O2 sheets as Na 
can be distributed over two sites that coordinate to one 
water molecule, while maintain an average coordination 
close to the ideal Na-6D2 0. 
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